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Abstract 

On the basis of the recently proposed Thermal Wave Model (TWM) for par- 
ticle beams, we give a description of the longitudinal charge particle dynamics in 
circular accelerating machines by taking into account both radiation damping and 
quantum excitation (stochastic effect), in presence of a RF potential well. The lon- 
gitudinal dynamics is governed by a 1-D Schrodinger-like equation for a complex 
wave function whose squared modulus gives the longitudinal bunch density profile. 
In this framework, the appropriate r.m.s. emittance scaling law, due to the damp- 
ing effect, is naturally recovered, and the asymptotic equilibrium condition for the 
bunch length, due to the competition between quantum excitation (QE) and radia- 
tion damping (RD), is found. This result opens the possibility to apply the TWM, 
already tested for protons, to electrons, for which QE and RD are very important. 
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In the study of charged particle beam dynamics for accelerators and plasma physics, 
a number of nonlinear and collective effects are relevant [0]. Due to the electromagnetic 
interactions between the particles and their image charges induced in the surroundings, 
these nonlinear effects also acquire collective nature This property is enhanced for 
very intense beams, which are employed in very high luminosity colliders. In addition, 
radiation damping and quantum electromagnetic fluctuations (quantum excitation) are 
generally present in the beam longitudinal dynamics, and, in particular, for electron 
circular accelerating machines are not negligible 0. 

Recently, a Thermal Wave Model (TWM) for charged particle beam dynamics has been 
formulated and successfully applied to a number of linear and nonlinear problems in 
beam physics In this approach, the beam transverse (longitudinal) dynamics is 

formulated in terms of a complex function, the so called beam wave function (BWF), 
whose squared modulus is proportional to the bunch density. This wave function satisfies 
a Schrodinger-like equation in which Planck's constant is substituted with the transverse 
(longitudinal) bunch emittance 0,^]. 

In particular, this model is capable of reproducing the main results of the conventional 
theory about transverse beam optics and dynamics (in linear and nonlinear devices) 0, 
and it represents a new approach for estimating the luminosity in particle accelerators 
0, 0, as well as to study the self-consistent beam-plasma interaction As far as the 
longitudinal bunch dynamics is concerned, in this scheme on can describe, in a simple 
way, the synchrotron motion when both self-interaction and the radio frequency (RF) 
potential well are taken into account. In particular, the right conditions for the coherent 
instability in circular machines have been recovered and soliton-like solutions for 

the beam density have been discovered 

In this letter we improve the thermal wave model for longitudinal bunch dynamics 
given in 0-0. By starting from the conventional longitudinal single-particle dynamics 
in circular accelerators, the problem under study is formulated in terms of an appropriate 
wave model which describes the evolution of the beam, when the RF potential well is 
taken into account together with radiation damping and quantum excitation. We show 
that the longitudinal beam dynamics is still correctly governed by a Schrodinger-like 
equation for the BWF. The envelope description is straightforwardly obtained from the 
wave solution and, correspondingly, the results are compared with those that are given 
in the conventional theory. In particular, an asymptotic time-limit for the bunch length 
and the r.m.s. emittance scaling law are obtained. 

Let us consider the motion of a single particle within a stationary bunch, travelling with 
velocity /3c (/? ^ 1) in a. circular accelerating machine of radius Rq = cTq/2'k {Tq being 
the revolution period). It is well-known that if both radiation damping and quantum 
excitation are taken into account, defining x as the longitudinal displacement of the 
particle with respect to the synchronous one, and s = ct > cTq {t being the time), the 
particle dynamics is governed by the following set of equations |TD[ : 
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where V = ^ denotes the relative longitudinal energy spread of the single particle, com- 
puted with respect to the synchronous one {AE = 0), and U {x, s) = (I/cTqEq) Jq q AV{y) dy 
is the effective potential energy that the particle sees after a turn in the ring {AV be- 
ing the corresponding total voltage variation seen by the particle). Moreover, in Eqs. 

Eq and q stand for the synchronous particle energy and charge respectively, and 
?7 = (1 — P^) — ttc is the phase slip factor is the momentum compaction) . 
Finally, in (|^), 7 represents the damping coefficient [|10| and dR/ds accounts for the 
quantum excitation effect (noise), R{s) being the difference between the energy effec- 
tively radiated by the particle during the time interval s/c and the average of this energy. 
Since R{s) is a stochastic quantity, namely its average value vanishes whilst its r.m.s. is 
not zero, we cannot treat the force term dR/ds as the other ones of r.h.s. in (Q), which on 
the contrary have deterministic nature. Hereafter, without loss of generality, we restrict 
our analysis to electron machines (g = e, \ri\ etc) for which the e.m. emission is a 
relevant phenomenon. 

The r.m.s. of R'^{s), R^, obeys to the following equation 

+ jR2 = , (3) 
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where po is the magnetic bending radius, < >^/^ is the quantum fluctuations associated 
to the noise, A^-^ is the mean rate of photon emission, rg is the classical electron radius, h 
is the Planck's constant, and 70 = (1 — /3^)~^/^ is the relativistic factor. 

It is easy to see that (|l]) and (Q) are the usual equations for the longitudinal motion 



with 



Ton under the substitution s = ct and V = AE/Eq. 

It is worth to point out that the distinction between deterministic and stochastic 
force terms is a necessary requirement to properly construct the wave equation ruling the 
collective behaviour of the beam. This difference can be understood observing how the 
emittance, which also describes a stochastic effect but related to the temperature of the 
system, is involved in a wave equation (Schrodinger-like equation) 1^-0. In fact, despite 
of deterministic terms, the emittance is the only quantity which plays the role analogous 
to a diffraction parameter. In the framework of TWM this parameter is involved in the 
quantization rules. 

To understand how to include quantum excitation in the TWM description, we first 
consider the simplified situation in which quantum excitation is negligible. Then, we 
generalize our results by taking into account also this effect. Under this hypothesis, by 
considering a linearized RF- voltage only {U{x, s) = Urf{x) ~ {K/2r]) x^, where K is the 
RF cavity strength, supposed for simplicity to be positive), it is easy to prove that the 
Lagrangian associated to (|l]) and (|]) is given by [|lT], |jl2l 



C(x, x', s) = — x'^ — K x^ e'^ , (5) 

where x' = dx/ds. Hence, computing the x- variable conjugate momentum hj p = 
dC/dx' = {x' /rf) exp(7s), the corresponding hamiltonian results to be 

n{x,p, s) = ^p^ e-^' + \—x^e^' . (6) 
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In order to write a Schrodinger-like equation for the BWF, which describes the longitudinal 
dynamics of a short bunch (a << Rq) in presence of the only radiation damping we have 
to follow the Thermal Quantization Rules (TQM), in complete analogy with our previous 
works [^-0 : p ^ p = —ied/dx , and H ^ H = ied/ds. Therefore, (|^) gives (for 

d v^e^ e"^''^ (9^ 1 
zr/e e"^^— *(a;,s) = — ^(x,s) + -Kx" ^!{x,s) , (7) 

where e is a constant to be determined, which, according to Ref.s accounts for 

thermal spreading of the bunch and plays the role fully similar to Planck's constant. The 
BWF \E'(x, s) satisfies the following normalization condition 

|^(x,s)p dx = I , (8) 

which, fixed for s = 0, holds for any s, due to hermiticity of the hamiltonian operator 
(^. Thus, according to the features of the TWM, if is the total number of particles 
in the bunch, N\^{x,s)\'^ represents the longitudinal bunch number density (number of 
particles per unity length). 

Interestingly, a complete set of solutions of (|^) is given in terms of Hermite-Gauss modes: 



{27r22™(m!)V2(s)}'/^ 1 4a2(s) j " \V2a{s) 
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where the are the Hermite polynomials (m = 0, 1, 2, ...), and the functions cr(s), p(s) 
and 0(s) satisfy the following set of differential equations 

(Pais) d ais) , , rfe^ e'^T^ , , 

+ -i—^ + Ka{s)- ' „. . =0 . (10) 



ds^ ' ds ^ ' 4 (t3(s) 

1 _ 1 da{s) 
p (j(s) ds 



. (12) 
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In particular, taking (P) for m = (fundamental mode), we note that cr(s) results to be 
the corresponding bunch length, defined as 



a'is] 



/ + 00 
\'^o{x,s)\'^ dx =< x'^ > . (13) 
-oo 



In addition it is useful to define the momentum spread ap{s 

d 
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dx =<p^ > . (14) 
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By using the definition of p and of \E'o(a;, s), and the quantum formahsm for the average 
of the operators, we can show that the following expression A{s) is a constant of motion 



Ms) 



< X > 



1 . . 
- < xp + px > 



1/2 



(15) 



and it coincides with the diffraction parameter e of the TWM. Observe that (|T5]) is 
formally identical to the well-known Robertson-Schrodinger uncertainty relation |T3|,[p!^, 
and e, which is straightforwardly obtained taking the minimum of this relation, is the 
natural extension to the quantum-like description of one of the Co urant- Snyder invariants, 
well-known in particle accelerators (Poincare-Cartan invariants, in classical mechanics 

El). 



From (|1^), we can also derive the scaling law for the following effective emittance e(s) 
{quantum-like r.m.s. emittance), defined in analogy to the classical definition of r.m.s. 
emittance given by LapostoUe [P]],|T7[| 

1/2 

(16) 
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The same scaling law can be extrapolated from Eq. (|TUp by comparing it with the cor- 
responding expression for the undamped case [§, consequently the diffraction parameter 
of TWM e represents the initial value e(0). Furthermore, by virtue of ([T6|) , the envelope 
equation (^) could be rewritten substituting eexp(— 7s) with the effective emittance e(s). 
It is worth to point out that the presence of the friction-like term in (p!0D (7 7^ 0) is re- 
sponsible for the s-dependence (time dependence) of r.m.s. emittance, whereas in the 
undamped case (7 = 0) e is a constant of motion and gives the accessible phase-space 
area associated with the system. 

Remarkably, in order to obtain all the above results concerning with radiation damping 
we could start from the following Schrodinger-like equation 



d 

ir]e{s) — ^(x,s) 



„2^2/ \ Q2 1 

^^£^^(x,.) + '-Kx'^{x,s) 



(17) 



formally obtained by introducing in (^) the following TQR: p ^ p = —ie{s)d/dx , 
and Ti. ^ Ti. = ie{s)d/ds, where e(s) is given by (|16|). For the present case we note 
that: 

a) e(s) satisfies the following differential equation 

de 

+ 76 



ds 







1^ 



with the initial condition e(0) = e; 

b) for this new TQR the Robertson-Schrodinger-like relation ([T5|) does not give a con- 
stant of motion anymore, since it results to be: A{s) = e e""^* = e{s). 



The above remarks allow us to immediately generalize our results to the case in which 
both deterministic (RF potential well plus radiation damping) and arbitrary stochastic 
effects are taken into account. In this case, in fact, the system dynamics is assumed to 
be ruled by 

nie{s) —m{x,s) = -^^^^(^'^) + ^Kx''m{x,s) , (19) 
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where now e(s) is an arbitrary function of s, to be specified in correspondence of the 
particular stochastic effects considered, but satisfying the initial condition e(0) = e. Also 
in this more general case of damped harmonic oscillator, we are able to give a complete 
set of normalized solutions of Eq. (|T9D 

I i-C I ^ ^ ( iAj 

'^m{x,s} = , ,,1/4 exp 



{27r22™(m!)2(T2(s)}'/^ 1 4cr2(s) j " VV2t^(s 



xexpit +z(l + 2m)0(g) , (20) 

I 2r7e(s) p(s) J 

where now the function a{s) satisfies the following equation 

with r(s) = — rflog[e(s)]/(is, and p{s) and 0(s) have the same definition of (|TT]) and ( [T^ ) 
in terms of o"(s). By defining a(s) as 

a(s) ^a(.) exp{^J\is')ds'] = ais) ^ , (22) 

Eq. (|2l|) becomes 



where ^ A' - - i^™. 

Some physical consideration are in order 

1) At the early time (7s << 1), e ~ e, and a{s) ~ 5"(s). We physically expect 
that, during this time scale, the damping rate r(s) is maximum whereas the quantum 
excitation is negligible. In fact, at the beginning, due to the small number of produced 
photons, the photon noise is negligible. Consequently, r(s) ~ 7 = const., and, thus, 
K[s) K — (7^/4) = const. In particular, for K > (7^/4), 7^ is a positive constant, so 
that a{s) and a{s) are limited functions. 

2) For very large time (7s >> 1), since in the usual description of the circular ac- 
celerating machines |^ a sort of asymptotic equilibrium is reachable, we can seek this 
equilibrium by assuming that r(s) vanishes as s increases its values. This physically 
means that, as s grows, the radiation emission produces more and more photons which, 
in turn, increase the quantum excitation. Since the radiation damping provides a decreas- 
ing of the r.m.s. emittance, a sort of competition between this damping and quantum 
excitation is thus established in such a way to reach an asymptotic limit for e(s), say 
ez). Consequently, from the asymptotic condition r(oo) = 0, it follows that K = K > 0. 
Thus, we conclude that, for 7s >> 1, a{s) is required to be limited, and from (P^ ) and 
(p^) it follows that asymptotic equilibrium solutions a{oo) and a{oo) for a{s) and <j{s), 
respectively, exist and results 



5-(oo) = (t(oo)J— . (24) 




Since the present description holds for an arbitrary form of e(s), on the basis of the 
above physical considerations we have to assume its explicit behaviour or, equivalently. 
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to assume the form of r(s), in such a way to recover the results of conventional theory 
1^. To this end, we observe that (|18D , which is valid in the case of negligible quantum 
excitation (e/? = 0), suggests to assume in the more general case eo 7^ 0, the following 
evolution equation for e{s) 

de ^ 

— + 7e = 7eD , (25) 
with e(0) = e. Consequently, we have 

e{s) = ee-^' + 6^3(1 - e"^^) , (26) 

and 

(e - eo) + eoe^" 
We observe that for 7s >> 1, ( p^D becomes 

trjen^^ix,s) = -"U^ ^m{x,s) + ]^Kx^^!{x,s) , (28) 

and the uncertainty relation at the minimum gives 

a(cx))ap(oo) = y , (29) 

where crp(oo) = lim^^oo < >• Note that o"p(oo) represents the expectation value of the 
longitudinal momentum spread of the bunch at the asymptotic equilibrium. This means 
that its explicit determination does not depend on the history which led the system to the 
equilibrium condition, but it depends on the quantum fluctuations (photon noise) only. 
Furthermore, (t(oo) can be explicitly determined by imposing the asymptotic equilibrium 



in the (0). Thus 

which, by using (|29|), becomes 



a\oo) = ^ , (30) 



a(oo) = ^— — ^ crp(oo) , (31) 

where z/^ (synchrotron number) stands for the ratio between Vt^ = c\fK (synchrotron 
frequency) and the revolution frequency cjq = Pc/Rq ~ c/Rq. Consequently, since from 
(H) the equilibrium value {R^)eq of is Qp/l, and since it is proved that the equilibrium 
energy spread is given by \{R^)eq 0, it is very easy to recognize that 



Hence, by combining (P^ and (^TJ) we finally obtain 

1/2 

(33) 
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Eq. (^) recovers the well-known proportionality between the equilibrium bunch length 
cr(oo) and the corresponding equilibrium momentum spread (Tp(oo), given in literature. 



7 



whilst (^2]) and (pSD give explicitly these values in terms of the quantum fluctuations. We 
remark that, as a{s) and <7p{s) go to the asymptotic equilibrium values, e goes to the 
minimum value eo of the r.m.s. emittance starting from the initial value e, according to 
the scaling law (|26|). We point out that the reduction of e(s) is due to the bunch cooling 
produced by the radiation damping, and, consequently, en represents the limit value of 
e{s) for which the bunch thermalization is completed. In addition, by combining (pQ]), 
(p^) and ( |33D we immediately get the expression of eo in terms of electron Compton's 
wavelength Ag = h/nieC {rrie being the electron rest mass) 



en 
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= XD 



(34) 



In table 1 we report the values of the factor xd for some electron-positron circular ma- 
chines. Note that the large values found for xd imply eo » K- This suggests that at 
the equilibrium the quantum-like behaviour of the system as a whole still corresponds to 
the Liouville regime (thermal equilibrium) and, according to (pPj), would represent a sort 
of macroscopical coherence. On the other hand, the quantum limit of such a coherence 
is, of course, recovered for e^) ~ Ag (Heisenberg regime). 



In this letter we have presented an extension of the recently proposed thermal wave 
model for particle dynamics ^ to the longitudinal motion in circular accelerating machines 
when both RD and QE are taken into account. In this framework, the particle dynamics 
in the presence of a RF potential well is governed by a 1-D Schrodinger-like equation for a 
complex wave function, whose squared modulus gives the longitudinal bunch profile. We 
have shown that the solutions for the BWF of this problem are given in terms of the well- 
known Gauss-Hermite modes. In particular, the fundamental mode (lowest-energy mode) 
gives a pure Gaussian space-distribution for the particles, and the corresponding envelope 
equation gives an asymptotic value for the bunch length, which is expressed in terms of 
the quantum fluctuations (noise). Correspondingly, according to (|^), as the beam is 
cooling, due to RD, the emittance goes to the equilibrium value eo which represents the 
equilibrium limit and plays the role of thermalization value. 

In conclusion, the above results allow us to apply the thermal wave model, already suc- 
cessfully applied to the undamped longitudinal dynamics (protons) [|| , to the synchrotron 
electron motion in a more accurate way, since in this case both radiation damping and 
quantum excitation are not negligible. This occurs, for example, for electrons in circular 
accelerating machines. 
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Table 1. 
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1 1 

1^1 ~ Oic 
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1.02 • 10 
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7.63 • 10 ^ 


3.74 • 10 ' 
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1.2 • 10^ 
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1.55 • 10-3 
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1.8 • 10^ 


2.44 • 10-=^ 


5.22 • 10-2 


3.6 • 10-^ 


3.4 • 10^ 


LEP 


3.1 • 10^ 


1.1 • 10^ 


2.9 • 10-^ 


8.9 • 10-2 


3.54 • 10-^ 
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